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The theory by Cross and Fisher (CF) is by now commonly accepted for the description of the spin- 
Peierls transition within an adiabatic approach. The dimerization susceptibility as the essential 
quantity, however, is approximated by means of a continuum description. Several important exper- 
imental observations can not be understood within this scope. Using density matrix renormalization 
group (DMRG) techniques we are able to treat the spin system exactly up to numerical inaccuracies. 
Thus we find the correct dependence of the equation of state on the spin-spin interaction constant 
J, still in an adiabatic approach. We focus on the pressure dependence of the critical temperature 
which is absent in the CF theory as the only energy scale with considerable pressure dependence is 
J which drops out completely. Comparing the theoretical findings to the experimentally measured 
pressure dependence of the spin-Peierls temperature we obtain information on the variation of the 
frustration parameter with pressure. Furthermore, the ratio of the spectral gap and the transition 
temperature is analyzed. 
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I. INTRODUCTION 

Low dimensional quantum systems are currently of 
considerable interest mainly due to the fascinating phase 
transitions driven by strong quantum fluctuations. The 
continuous interest from the theoretical side is provoked 
by the discovery of many experimental systems realizing 
quasi one-dimensional quantum systems. In the field of 
spin-Peierls systems the discovery of the inorganic com- 
pound CuGeOs realized a milestone as many measure- 
ments have been performed with high accuracy since. 
Therefore, CuGeOa has attracted much attention in ex- 
perimental as well as in theoretical works. The high 
temperature behaviour of CuGeOa was found to be mod- 
elled adequately by one-dimensional frustrated Heisen- 
berg chains In the dimerized phase, many features 

were shown to be consistent within an adiabatic descrip- 
tion of the phonon degrees of freedom. This observation 
comprises zero temperature |5|-]^] as well as thermody- 
namic properties 

Even in one space dimension only a few exact results 
exist particularly concerning thermodynamics. For in- 
tegrable systems the thermodynamical potentials and 
asymptotic behaviour of correlation functions are known. 
A notorious problem is posed by response functions and 
non-integrable systems in general. With respect to this, 
the recently developed transfer matrix DMRG (TMRG) 
[p[-|ri| on the basis of transfer matrices provides a 
very powerful method to calculate thermodynamic quan- 
tities of spin chains without any use of perturbative meth- 
ods. This has been demonstrated in several applications 

In this paper we study the influence of microscopic 
coupling constants on the spin-Peierls transition tem- 
perature. This allows for an understanding of the 
considerable pressure dependence of the phase diagram 



along the following line of reasoning. It is known 
that external pressure affects the magnetic properties 
of CuGeOa considerably |^,^,^. Fits for the magnetic 
susceptiblity yield the change of the nearest-neighbour 
spin interaction J and estimates for the frustration 
paramter a. Using these data for J as function of pres- 
sure we are able to explain the observed increase of the 
spin-Peierls temperature and estimate the pressure de- 
pendence of the next-nearest-neighbour exchange 

The outline of the paper is as follows. In section |l| we 
present the model and a motivation for our description 
of the experimentally studied spin-Peierls systems. We 
study the static dimerization susceptibility in section III . 
Section IV is devoted to the computation of the critical 
temperature as function of the spin exchange couplings 
and external pressure, respectively. We give a compari- 
son of our results with experimental measurements. In 
section ^ we investigate the spectral gap and its ratio to 
the spin-Peierls temperature. The conclusion is given in 
section 

II. MODEL 

In the inorganic spin-Peierls compound CuGeOs the 
magnetic interactions are attributed to Heisenberg spin 
exchange. There is numerous evidence that in addition 
to the nearest-neighbour interaction (J) a next-nearest- 
neighbour exchange J' — aJ with a — 0.35 has to 
be taken into account. Usually the constant a is referred 
to as frustration parameter. 

At the spin-Peierls transition the system undergoes a 
structural phase transition driven by the quantum spin 
system coupled to the phonons. The spin-phonon cou- 
pling is modelled by spin exchange integrals depending 
linearily on the local displacements. The adiabatic treat- 
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ment yields a quantum spin system coupled to just one 
phonon mode with the commonly used Hamiltonian 



H 



E 



J[(1 + ,50S,S, 



K 



(1) 



where S,; are spin 1/2 operators, 5i ~ (— denotes 
the modulation of the magnetic exchange couplings in 
the dimerized phase. Here, we restrict ourselves to van- 
ishing external magnetic field where the system shows a 
phase transition from the uniform (U), i.e. (5 = to the 
dimerized (D) phase ((5 > 0). 

The elastic energy can be expressed in terms of mi- 
croscopic constants rendering equivalent to an RPA 
treatment of the phonon propagator for the full spin- 
phonon system. Within RPA the condition for the phase 
transition is identical to that for (|l|) as formulated in (Q) 
below if K is adjusted to the following value pffl 



K 



A 



(2) 



The sum runs over the spin-Peierls active modes. For 
each mode i, denotes the effective mass of the unit 
cell, the frequency of the spin-Peierls active phonon 
and Xi the spin-phonon coupling constant. In particular, 
K is proportional to J^. The constant C contains only 
the microscopic parameters of the underlying lattice. 

Our numerical investigations do not improve over those 
of CP with respect to the RPA treatment. However, 
within the RPA approximation we deal with the complete 
dynamics of the quantum spin system. Note that in Ref. 
pl[ a continuum description of the spin system was used 
with a subsequent bosonization treatment which is be- 
lieved to capture only the long distance asymptotics of 
the correlation functions. 



III. RESPONSE FUNCTIONS 

The static dimerization susceptibility of the spin sys- 
tem is defined by 



AJx) 
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where x = T/ J, and fa{x, 6) is the free energy per site for 
system with fixed dimerization i5, frustration a and 
K set to zero. The response function A is nothing but 
the correlation of the nearest neighbour spin exchange 
SiSi+i (dimer operator) at momentum q = ir and energy 
uj — 0. The U/D phase transition takes place for 



Aa,{xsp)^K/J = CJ. 



For details the reader is referred to 



(4) 



Let us now review the results obtained by CF. Within 
the bosonization approach they find x ■ Aa=o{x) = xo 
with Xo ~ 0.26. As a direct consequence by use of the 
inversion of (Q), 



Tsp = JA-\CJ), 



(5) 



this yields dT^p/dJ — 0. Of course, this is also clear 
from the fact that the energy scale of the spin system 
completely drops out due to scale invariance. 

Fig. shows a comparison of our TMRG rcsultf^ for 
the function x ■ Aa{x) for various values of a with the 
findings of CF and exact data for free fermions. The 
enormous progress achieved by the numerical analysis 
is the correct treatment of the spin system on the lat- 
tice at practically all length scales. This improves over 
the continuum limit approach in which the asymptotics 
of correlation functions is incorrectly extended to short 
distances. For the unfrustrated Heisenberg model, i.e. 
a = 0, we are able to observe directly the deviations be- 
tween the continuum limit and a lattice treatment of the 
spin degrees of freedom. 



xA(x: 




FIG. 1. Depiction of the function x-A{x) with x — T/J for: 
free fermions, Heisenberg model in continuum limit (CF), and 
TMRG results for frustration parameter a = 0,0.241,0.35. 
The circels denote the relevant values for CuGeOs (see text). 

With respect to CuGeOs we fix J by the requirement 
that the experimental magnetic susceptibility equals that 
of the strictly one-dimensional model at the critical point 
leading to J = 130K (a = 0), J = 150^ (a w 
ttc ~ 0.2412 0j2|,|2|), J = 1607^ (a = 0.35), and 
J — 350K (free fermions). The circles in Fig. [| denote 



^We have used 24 states in the renormalization step and 
have set the accuracy with respect to Trotter decomposition 
to (rM)-i = 0.05. 
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the values of ) at the experimentally determined 

spin-Peierls temperature Tgp = lAAK for CuGeOa (see 
e.g. p^-p6t). These values imply constants (of. Eq. ||) 
K = 0), Co.24i « 0.040 (a = 0.241), 



(a ^ 0.35), and C„. « 0.0027 K"' (free 



Co « 0.019 K"' 
Co.35 « 0.070 K- 
fermions). 

For the unfrustrated model the value xsp ■ ^Q=o(a;sp) 
almost coincides with the results of CF. However the 
agreement happens only fortuitously since xo is a zero 
temperature quantity. For other frustration parameters 
qualitative and quantitative deviations from the CF-line 
appear. The divergence oi A{x) for free fermions, Heisen- 
berg model with a < a — ac, and a > Ofc, is log(a;). 



l/a;xlog. corrections 
therein), and exponential 




1/x (see and references 
, respectively. The quanti- 



tative values for xa,p • A{xsp) also differ from each other. 
We must conclude that it is risky to deduce quantita- 
tive results from the bosonization approach as already 
pointed out by CF. For applications to CuGeOait is fur- 
thermore uncertain if the pure Heisenberg chain can even 
yield the qualitative results correctly. From the analysis 
of the susceptibility data at higher temperature we are 
led to favour the frustration parameter a = 0.35 |^-^] 
for which the behaviour of the response function devi- 
ates considerably from x^jx. 



IV. PRESSURE DEPENDENCE OF THE 
MAGNETIC SYSTEM 

From (H) it is straight forward to deduce the depen- 
dence of the spin-Peierls temperature on the variation of 
the spin coupling constants. The relation is valid for ev- 
ery fixed value of a. We first focus on the system with 
ao = 0.35 at ambient pressure. 

The dependence of J and a on the external hydrostatic 
pressure have already been obtained from the relation 
between magnetostriction and the pressure dependence 
of the magnetic susceptibility x ^J. From our TMRG 
data we find ^ d\aa constant temperature. 

These two quantities appear in the following expression, 
f = afe^, where x^ ^ J and x^ = a. There- 
fore, we conclude the value for the pressure dependence 
of J to be more reliable than that of a. The authors of 
Ref. @ deduced a value of ^ = -7.0(5)%/GPa us- 
ing the relation between magnetostriction and pressure 
dependencies of J and a. Using in addition the exper- 
imentally determined value for the pressure dependence 
of the spin-Peierls temperature « 4.8 K/GPa fTl] 
we obtain the relation T'sp(>^) (thick black line in Fig. 0). 

Due to the particular geometry involved in the super- 
exchange mechanism we expect the magnetic exchange 
energies to respond much more sensitively to the pres- 
sure (bond-bending mechanism p8| , p9| ) than the phonon 
frequencies or the spin-phonon coupling constants. We 
therefore consider C as independent of pressure which 



applicability for CuGeOa will break down at higher pres- 
sure. The numerical results with C = Co.35 are shown in 
Fig. ||. Obviously, a constant value of a = 0.35 (dashed- 
dotted line) can not explain the observed behaviour, even 
yielding a change of the critical temperature to opposite 
direction. 
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FIG. 2. Variation of the spin-Peierls temperature versus 
the relative change of the magnetic exchange coupling. The 
lines show TMRG results for constant a and fixed Co.35. The 
thick black line displays the experimental behaviour as an 
implicit function of pressure. The dashed (dashed-dotted) line 
shows the theoretical results for dj' /dp = {da /dp = 0). The 
corresponding dependencies of a on the pressure are shown 
in the inset. 

From geometrical reasons we expect J' to be non- 
decreasing. Assuming J' = constant to be realized, the 
increase of Tsp is too small by a factor of approximately 
4 (as displayed by the thick dashed line in Fig. ^). Con- 
sequentially, the main effect must be a strong increase 
of J'. The deduced dependence of a on the hydro- 
static pressure fits well with = 24 ± 2%/GPa (or 

d^-^J' ^ 17 ± 2%/GPa respectively) up to about 1.4 GPa 



dp 



as displayed in the inset of Fig. 
mined assuming that the value of 

±0.5-fi'/GPa, where the uncertainty in 



. The error was deter- 
involves an error of 

In J 



dp 



only playes a 



secondary role. In general, already the weak requirement 
that J is non-increasing with hydrostatic pressure gives 
a minimum pressure dependence of > 20%/GPa. 



The pressure dependence of a has already been investi- 
gated in Ref. 1^, however, with large error bounds which 
are respected by our results. The divergence of 
at a finite pressure or in other words, an upper limit of 
the critical temperature, is a physical prediction of the 
chosen one-dimensional approach. It is mostly based on 
a limited spontaneous gap as a function of a. But as 
already mentioned above, we expect an agreement with 
CuGe03 only in the low pressure region. 
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A similiar analysis can be done on the basis of the 
unfrustrated Heisenberg model and free fermions. Eval- 
uating the magnetic susceptibility data of Takahashi et 
al. @ we derive ^ « -5%/GPa (a = 0) and 

« — 6%/GPa (free fermions). In contrast to the 
frustrated change of a under pressure for these 

initially unfrustrated models is not reasonable. The de- 
duced theoretical BCS-ratio is too small by a factor of 
about 3 for free fermions (or equivalently Hartree-Fock 
calculations), in the case of the unfrustrated Heisenberg 
model even by a factor of « 15 (see Tab. I). 

We must conclude that neither the unfrustrated 
Heisenberg model nor Hartree-Fock results are able to 
describe the physics of CuGeOs correctly. In constrast 
to this a Heisenberg model with parameters uq = 0.35 
24%/GPa turns out to reproduce the exper- 



dp 



imental findings up to about 1.4 GPa. 



V. "BCS-RATIO" 

Another interesting quantity is the ratio of the singlet- 
triplet gap to the spin-Peierls temperature, known as 
"BCS-ratio" . Combining the TMRG data for A with zero 
temperature DRMG calculations for the singlet-triplet 
gap as a function of a and the dimerization one gets a 
relation between Ag^ and Tsp- The results are shown in 
the bottom section of Fig. ^. 
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FIG. 3. Bottom: Singlet-triplet gap Ast as a function of 
the spin-Peierls temperature. Symbols denote the DMRG re- 
sults, the solid line shows the exact result for free fermions. 
The thin lines show the expected behaviour from the scaling of 
A (see text). Top left: Ratio of the spin gap and spin-Peierls 
temperature (BCS ratio). Top right: BCS ratio, experiment 
versus theory. 

Using the scaling of ^ the definition of the criti- 
cal temperature and the dependence of the ground state 



energy on small saturation dimerizations Sq = 5{T = 0) 
one finds 



(6) 



ignoring logarithmic corrections for a < ac. Here Ag.^, de- 
notes the singlet-triplet gap for vanishing dimerization, 
which is zero for a < ac- In the next step we apply 
Ast — Ag.j, ~ 6q^^ ^^^^ f o derive an explicit asymptotic 
relation between the gap and the spin-Peierls tempera- 
ture, 

AO 

Ast - AO, c Tsp exp(-^/rsp), (7) 

again neglecting logarithmic corrections for a < ac- The 
thin lines in Fig. || show fits according to Eq. (Q) in the 
range up to T/ J = 0.15. For a = 0.5 we used the known 
value Ag.j, = 0.2338 J |^ , ^ , hence only performed a 
one-parameter fit in the presented region. For a = 0.35 
a spontanous gap of A^.^, = 0.035 J is used, also derived 
by T = DMRG. For a = the predicted hnear be- 
haviour (cf. (0)) can not be seen due to the logarithmic 
corrections. A linear extrapolation of our data points 
shows a positive offset. However, this is consistent with 
the presence of logarithmic corrections since we expect a 
zero limit at T = with infinite slope. The free fermions 
show the well known behaviour Ast/Tsp ~ 1-76. We now 
derive the BCS-ratio as a function of the spin-Peierls tem- 
perature simply by dividing by Tsp (upper left of Fig. |^). 
Fistly we observe that the Heisenberg like models have a 
distinctly larger BCS-ratio than free fermions. 

For comparison to experiment we again fix the con- 
stant C. Furthermore, the pressure dependence of J and 
the frustration derived in section IV is taken into ac- 
count here. The results are shown in the top right of 
Fig. 1^. Interestingly, the experimental result ||l8|| com- 
pares well with the low temperature asymptotics for free 
fermions but showing a pronounced increase with pres- 
sure which is not reproduced. From the investigations in 
Ref. [H we already know for Heisenberg models that the 
gap and therefore the BCS-ratio is larger than the exper- 
imentally observed one. This happens due to the strictly 
one-dimensional treatment of CuGeOs , i.e. the neglect 
of the dispersion perpendicular to the chain which will 
lower the true gap. The unfrustrated chain even yields 
a qualitatively incorrect tendency of a decreasing BCS- 
ratio under pressure the appropriately frustrated system 
is at least able to explain the increase. 



VI. CONCLUSION 

The TMRG analysis of the spin-Peierls phase transi- 
tion allows a complete treatment of the quantum dynam- 
ics. In contrast to continuum descriptions correlations 
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are respected at all length scales, which leads to the ex- 
act response functions. Within the scope of an adiabatic 
description, equivalent to RPA as in the CF theory, we 
are able to study the influence of the spin-spin exchange 
energy scale J on the critical temperature. We like to 
emphasize that these results are a non-trivial improve- 
ment over the CF theory which shows no dependence on 
J at all. Moreover, frustrated models can be investigated 
as well. 

Using the pressure dependence of J it is possible to 
study the dependence of the spin-Peierls temperature on 
pressure. Neither Hartree-Fock calculations nor the un- 
frustrated Heisenberg chain yield the strong increase as 
measured. Once again, our investigations favour a frus- 
tration of a = 0.35 for CuGeOa at ambient pressure. 
We find a rather strong dependence of the frustration 
on pressure, = 24 ± 2%/GPa, which agrees with 



earlier studies 

The analysis of the "BCS-ratio" also gives a clear in- 
dication that frustration is present in CuGeOa , even 
though some quantitative deviations can only be ex- 
plained by residual perpendicular couplings. There is 
evidence for a strong dependence of a on pressure from 
the "BCS-ratio" . 
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^[%/GPa] 


^[%/GPa] 


^[K/GPa] 


free fermions 


^ -6 


(0) 


^1.8 


Q = 


^ -5 


(0) 


^0.3 


ao = 0.35 


-7.0 ±0.4 


24±2 


4.8 ±0.5 



TABLE L Pressure dependences of the spin system param- 
eters and the spin-Peierls temperature as explained in the 
text. Numbers in bold face denote results obtained by our 
DMRG calculations, all other quantities were used as input 
data. Note that models based on free fermions or unfrustrated 
Heisenberg chains can not explain the rather large pressure 
dependence of the critical temperature. 
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